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Abstract 

We consider some of the issues surrounding the formation and evolution of Q- 
balls in the MSSM and its extensions. The ratio of the baryon number packed 
into Q-balls to that outside, fg, plays a fundamental role in determining the 
relationship of dark matter to the baryon number in the Universe at present. The 
final value of fs will depend upon the details of the formation of Q-balls from 
the collapse of the Affleck-Dine condensate and upon the subsequent evolution of 
the ensemble of Q-balls. We discuss the implications for neutralino dark matter 
in the gravity-mediated scenario and show that a light neutralino is necessary 
in most cases to account for the baryon to dark matter ratio, with an NMSSM 
singlino LSP of mass m x < 20 GeV being a favoured candidate. 



1 mcdonald@physics .gla.ac.uk 



1 Introduction 



The minimal supersymmetric (SUSY) standard model (MSSM) |l[] and its extensions 
offer a number of possibilities for the origin of the baryon asymmetry. In the MSSM, 
the most natural possibilities are electroweak baryogenesis || and Affleck-Dine (AD) 
baryogenesis ||. If we include the possibility of massive Majorana neutrinos we can 
also have leptogenesis [[|. 

AD baryogenesis rather generally results in the formation of Q-balls of baryon 
number (also known as B-balls) in the early Universe [5-11]. This occurs because 
the scalar potential of the squarks generally has attractive interaction terms [ O , 
making the homogeneous Bose condensate of squarks unstable with respect to spatial 
perturbations. As a result, the condensate fragments to eventually form Q-balls. 

In the conventional gravity-mediated MSSM, AD baryogenesis may be character- 
ized by the dimension d of the non-renormalizable terms responsible for lifting the flat 



directions [13[]. For R-parity (R p ) conserving models (necessary to eliminate renormal- 
izable B and L violating operators the dimension must be even. It has recently 
been shown that the d = 4 models are generally ruled out by the effect of thermal 
corrections to the AD potential Ul~4 |, which suppresses the baryon number. Thus d = 6 
models are favoured. For CP violating phases 5c p (expected to be of the order of 1) 
these must have a reheating temperature Tr « 1 GeV/5cp in order to account for the 
observed baryon asymmetry ||16|| . (Such low reheating temperatures can be a natu- 
ral feature of SUSY inflation models 0].) The resulting Q-balls are very long lived, 



typically to temperatures 1 MeV - 1 GeV JO]]. As a result, it is possible that the late 
decay of Q-balls to baryons and neutralinos can explain the remarkable similarity of 
the number densities of baryons and dark matter particles for the case of WIMP dark 
matter with mass 0(mw), fc> r example neutralinos |K], [TJJ. For Q — 1 and Qdm ~ 0.4, 
the ratio of the number density of baryons to the number density of dark matter par- 
ticles is given by n B /n DM « (1.5 — 7.3) m DM /m w [jT^|. (In the following we refer 
to this as the baryon to dark matter ratio.) This is a major motivation for the idea 
of late-decaying Q-balls as a simultaneous source of the baryon asymmetry and the 
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dark matter density. More generally, late-decaying Q-balls will provide a source of 
non-thermal WIMP dark matter, resulting in different predictions for the dark matter 
density as a function of the parameters of the MSSM HIT] , HH- This variant of AD 
baryogenesis, with late-decaying Q-balls carrying baryon number, is known as B-Ball 



Baryogenesis [11 



In the case of gauge- mediated SUSY breaking [5-9], where Q-ball formation from 
the AD condensate was first proposed [§, |6], [7[ , the potential of the AD scalar becomes 
essentially completely flat for scalar expectation values larger than the mass of the 
messenger fields which transmit SUSY breaking ||. As a result, the mass of the 
AD scalars in the Q-balls satisfies m oc Q -1 ^ 4 , and for large enough Q-ball charge 
the Q-balls are completely stable with respect to decay to nucleons and may account 
for dark matter |7|, [19|, with interesting experimental and astrophysical consequences 
[§]. The cosmology of AD condensate fragmentation in models of gauge-mediated 
SUSY breaking will be significantly different from the case of gravity-mediated SUSY 
breaking, but we expect aspects of the gravity-mediated case to be relevant to the 
gauge-mediated scenario. In this paper we will focus on the more conventional gravity- 
mediated SUSY breaking scenario. 

In both the gravity- and gauge-mediated scenarios the final baryon to dark matter 
ratio is crucially dependent upon the ratio of baryon number trapped in the Q-balls to 
that outside, Both scenarios require that /# < 1. Experimental limits from LEP 
and Tevatron require that the MSSM lightest SUSY partner (LSP) neutralino mass, 
m x , satisfies m x > 46 GeV POfl . (This assumes universal A-terms and gaugino masses 
at the unification scale.) As we will show, this implies that fs ~Z 0.4(fi x /0.4) when the 
Q-balls decay. In the gauge-mediated scenario, /b = 1 would result in all the baryon 
number today being in the form of stable Q-balls, an obviously unacceptable scenario. 

There have been some attempts to estimate fs by studying the classical dynam- 
ics of the collapse of the AD condensate to non-topological solitions. Enqvist and 
McDonald |21| considered the collapse of a single spherically symmetric condensate 
lump. It was observed that the value of fs depended upon the charge density of the 
original AD condensate relative to the maximum possible charge, which we denote by 
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Q. As we will show, if Q is less than one then not all the energy in the original AD 
condensate can be accomodated in the form of positively charged Q-balls. One possi- 
bility is that the condensate collapses to higher energy objects, Q-axitons [2lJ, which 
are classically stable but which can in principle decay to eventually reach a Q-ball. 
For Q-axitons it was found that Jb could be significantly less than one if Q <C 1, for 
example Q = 0.01 resulted in fg as low as 0.3 pi). On the other hand, in a lattice 
simulation of condensate fragmentation, Kasuya and Kawasaki |23| did not find evi- 
dence for fs significantly less than one. In addition, they also pointed out that it was 
possible to form positive and negative charged Q-balls, which could in principle allow 
all the energy in a Q < 1 condensate to be entirely accomodated in Q-balls. 

In this paper we will discuss various aspects of the /b problem in more detail. 
The question of whether ±Q-balls or Q-axitons form will be seen to depend upon the 
perturbation of the phase of the AD condensate; if the perturbation of the phase of 
the condensate is larger than the average phase when the condensate fragments, then 
±Q-balls will tend to form; if not, then Q-axitons will form. In the case where Q- 
axitons form there is the queston of how the classically quasi-stable Q-axitons evolve 
to lower energy Q-balls. We will suggest that this might occur classically via a slow 
emission of energy in scalar field waves or via annihilations of scalars within the Q- 
axiton. If the excess energy in the Q-axiton relative to the Q-ball is radiated before 
the Q-axiton decays and before thermal neutralinos freeze-out then it is possible that 
fs can be small enough to account for the baryon to dark matter ratio for a wide 
range of neutralino masses. 

A very different picture emerges if the Q < 1 AD condensate collapses to ±Q- 
balls. In this case, when the Q-balls decay, assuming /b ~ 1 for Q-balls as implied by 
numerical simulations, the decaying ±Q-balls will produce oppositely charged baryon 
number but contribute equally to the density of neutralinos, resulting in a large number 
of neutralinos relative to baryons for Q C 1. (A similar situation arises if Q-axitons do 
not lose their excess energy before they decay.) We will show that it is then essential 
that the MSSM neutralinos annihilate with each other to give an acceptable density. 
As a result, there will be no direct relationship between the number of baryons and of 
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dark matter particles. However, we will still have a non-thermal relic density of MSSM 
neutralinos from Q-ball decay, which is an interesting prediction in itself. (There may 
be some evidence for non-thermal dark matter from non-singular galactic halos |24| . ) 

If we have a maximally charged AD condensate, Q = 1, then there will be direct 
formation of Q-balls without a Q-axiton stage. In this case fg ~ 1 and MSSM 
neutralinos consistent with experimental constraints cannot account for dark matter 
directly via Q-ball decay. Annihilations must reduce the number of neutralinos, again 
losing the direct connection with the baryon number density although still producing 
non-thermal neutralino dark matter. 

Thus it will be seen that it is generally difficult for the baryon to dark matter ratio 
to be explained via Q-ball decay in the context of the MSSM, with the only possible 
exception being the limiting case where Q-axitons form and evolve to Q-balls before 
neutralinos freeze out. However, extensions of the MSSM with weaker lower bounds 
on the LSP mass, in particular the next-to-minimal SUSY standard model (NMSSM) 



25 j with light singlino LSPs p6| |27], , will be seen to be consistent with the baryon 



to dark matter ratio. 

The paper is organized as follows. In Section 2 we consider the classical value of 
fs from condensate fragmentation. We consider the initial conditions for the linear 
evolution of the perturbations and the conditions for the formation of Q-axitons versus 
±Q-balls. We then comment on the numerical simulations of condensate fragmentation 
and point out that a sufficiently long evolution will be necessary to reach the time when 
absorption/re-emission of scalar field waves from collapsing condensate lumps becomes 
negligible and the condensate lump can reach their equilibrium states, allowing Jb to 
be extracted. In Section 3 we discuss the possible evolution of a Q-axiton to a stable 
Q-ball. In Section 4 we consider the consequences for neutralino dark matter in the 
MSSM and the NMSSM. In Section 5 we comment on the effects of more complicated 
aspects of realistic condensate fragmentation on the value of fs- In Section 6 we 
present our conclusions. 
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2 Classical Evolution of the AD Condensate 



In this section we consider some aspects of the classical evolution of perturbations. 
The formation of Q-axitons versus ±Q-balls will depend upon the initial size of the 
perturbation in the phase of the AD condensate when coherent oscillations begin and 
the baryon asymmetry is established at H w m, where m is the mass of the AD scalars. 

2.1 Energy in Q-balls for a non- maximally charged conden- 
sate 

We first show that for a non-maximally charged condensate in the gravity-mediated 
scenario it is not possible for all the energy in the condensate to be accomodated in 
positively charged Q-balls. The energy density in a coherently oscillating scalar field 
condensate with amplitude <p is p = m 2 (f) 2 . The maximum possible charge density 
in the condensate is pQmax = p/m = mcpl, corresponding to the case where all the 
scalars in the condensate are charged. The maximum energy density in Q-balls is 
Pq-baiis = m Pq = m Q Pq max = Qp, since the mass of the scalars in the Q-balls is 
to a good approximation m in the gravity-mediated scenario [O and we assume that 
all the charge is packed into Q-balls. So the energy density in Q-balls to the energy 
density in the condensate is 

P^i = Q. (1) 
P 

Thus for a non-maximally charged condensate with Q < 1 not all the energy density 
can be accomodated in positively charged Q-balls. One must either form higher energy 
non-topological solitons (Q-axitons) or have an ensemble of positively and negatively 
charged Q-balls. 



2.2 Conditions for Q-axiton versus ±Q-ball formation 



We next review the evolution of perturbations at H > m [pi . Perturbations of the 
AD condensate may be expected to arise from quantum fluctuations of the AD scalar 
during inflation, with 5(f) ~ H/2n. The subsequent evolution of the perturbations will 
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depend upon the paramaters of the scalar potential. In general, the perturbations 
satisfy an equation of the form |21| 



+ 3H5<p - V 2 (50 = -kH 2 5ct) , (2) 



where k is determined by the parameters c, \a\\ and d of the AD potential pl| , 
Here K < is due to radiative corrections (with \K\ ~ 0.01 — 0.1) |10, 11, 1^], d is 



the dimension of the non-renormalizable term in the superpotential which lifts the flat 
direction (with d = 6 favoured), cH 2 gives the order H 2 correction to the scalar mass 
[p9H (with c positive and typically of the order of one for AD scalars) and we assume 
that the natural scale of the non-renormalizable terms is the supergravity mass scale 
M* = Mpi/\^8n 0. The A-term also receives order H corrections, A\ = A\ Q + a\H, 
where A\ is the gravity-mediated soft SUSY breaking term and a\ depends on the 
nature of the inflation model; for F-term inflation \a\\ is typically of the order of one 



H| whilst for minimal D-term inflation models it is zero 

For the case of the AD scalar, $ = (0! + i(j) 2 )/\/2, we find different values for k, 



ki and k%, corresponding to the equations of motion for 5(pi and 502 ||21|| . Thus the 
evolution of the S4>i and S(p2 perturbations can be quite different. As a result, the initial 
501 and <50 2 at if ~ m depends upon the parameters of the potential. In particular, 
if \a\\ is larger than |c|, then ki can be large compared with k\, which means that 502 
can be much more suppressed by expansion than 5<j)\ |[21|| . (During matter domination 



50 oc a~ Re( - v \ where r] = (3/2 - ^9/4-4jfe)/2, whilst during inflation <50 oc a~ Re{<7 \ 
where a = (3 - V9-4A;)/2 |2lJ.) Thus depending on |c| and \a\\, <50i ^> <50 2 or 
(50i -C (502 at H pa m are possible. 

As H becomes smaller than m, the A-term becomes dominated by A\ . The 
real direction is determined by the A-term, so if we define 0i and 02 relative to the 
real direction then these fields are rotated at H < m relative to their directions at 
H > m. (At H > m, <p2 is damped to zero due to the A-term in the potential pl| .) In 
addition, due to the A-term, the coherent oscillations of 0i and 02 will have a relative 



phase shift typically of the order of 1, resulting in a baryon asymmetry approximately 
proportional to the phase 9 = 02 /4>i o of the condensate field, where 0i and 02 o are 
the amplitudes of the coherent oscillations. If the 9 between is small compared with 
one and 50i 3> 502 at H ^ m, then 502 at if ~ m after this rotation will be given 
by 502 ~ #50i. This is the smallest initial value of the 502 relative to 50i possible at 
H w m. This limiting case is relevant to the possibility of forming Q-axitons rather 
than ±Q-balls and we will focus on it in the following. 

We can then consider the growth of perturbation until they become non-linear. 
This has been studied analytically for the case of a maximally charged (MAX) Q = 1 
AD condensate |7], [11], [Hp. However, it has not been considered in detail for the case 
of a non-MAX condensate, for which the amplitude of 0i is large compared with that 
of 02. The equations of motion are 



& + 3if0, - V 2 0* = -m 2 (l + K)<Pi - Km 2 Uog \ 2 , i = 1, 2 . (4) 



For 0! (t) 3> 02 (t) during most of the oscillation cycle, the equation for 0i will be similar 
to the case of the MAX condensate, except that 2 + 0| will be varying periodically 
in time rather than constant. Thus the equation for the growth of space-dependent 
perturbations of 0i will be similar to the MAX condensate and we expect that the 
0! perturbations will go non-linear (50i/0i > 1) essentially as in the case of the 
MAX condensate |7], [IT| . The equation for perturbations in 2 will, however, be 
different, since the log term is dominated by 2 . Even if 502 > 02 {t), the equation for 
perturbations of 0i(= 0i(t) + 50i) is not significantly altered so long as 50 2 /0i < 1. 
Thus the perturbation in the energy density (which is essentially determined by 0i) 
will not go non-linear until 50i/0i ^ 1, even if 502/02 ^ 1- 

We next show that 502/01 ^ 50i/0i at all times if 502 ^ 50i initially at H w m. 
To prove this we must show that the rate of growth of 502 is no greater than that of 
50i. We consider 0i ^> 2 , as is true for most of the oscillation cycle, and expand the 
equations of motion, 




50i + 3F50i - V 2 50 : 



m 2 (l + 3ir)50i - 2ifm 2 02502/0i , 



(5) 



502 + 3#50 2 - V 2 50 2 



m 2 (l + 7T)502 - 2fsrm 2 2 50i/0 



(6) 
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Since 95 (pi = (0 2 /0i)50i 1z 502, the magnitude of the \K\ dependent term of the 50 2 
equation, which is responsible for the growth of the 50 2 perturbations, is less than that 
of the 501 equation. Thus if 50 2 ^ 50i initially at H m m then at all times 502 ^ 50i. 
In particular, 50 2 ~ #50i is a consistent soultion of the equations. Therefore for 
502 "Z 50i at the onset of coherent oscillations, non-linearity will occur only once 50i 
goes non-linear, at which point the condensate will begin to fragment to condensate 
lumps. 

In general this means that 50 2 /0 2 ^> 9 is possible when the energy density in con- 
densate goes non-linear (50i ss 0i) and the condensate fragments i.e. the fluctuation 
in the phase of the condensate can be large compared with the average phase. In this 
case we expect ±Q-balls to form, since the charge density of the non-MAX condensate 
is proportional to 0i0 2 ~ 0i50 2 and 50 2 is periodic in space with positive and negative 
values. However, in the limiting case where 50 2 m 95<pi at H « m (and 9 is small 
compared with 1), 50 2 ps 9 5 (pi is a consistent solution throughout. Therefore when 
the condensate fragments we will have 50 2 /0 2 ~ 1 i.e. the fluctuation of the phase 
of the condensate is approximately the same as the average phase. In this limiting 
case we expect that Q-axiton formation will be favoured over ±Q-ball formation. A 
full numerical evolution of the perturbations will be necessary to determine exactly 
what happens in this case, but we expect that there will be dominant formation of 
Q-axitons together with a few ±Q-balls. 

So the question of whether Q-axitons or ±Q-balls form will be detemined by the 
fluctuation in the phase at the onset of condensate fragmentation. This will all depend 
upon the evolution of the perturbations after they leave the horizon during inflation, 
which in turn depends on the parameters of the AD scalar potential. We believe 
that both cases are possible, with ±Q-ball formation being typical but with Q-axiton 
formation as a limiting case. 
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2.3 Numerical simulations and absorption/re-emission of scalar 
waves 

We next consider the classical evolution of the non-linear condensate lumps. For the 
case of Q-balls (MAX condensate or non-MAX with ±Q-balls), this is straightfor- 
ward. Non-relativistic Q-balls carrying essentially all the charge form directly from 
the condensate and move apart from each other due to the expansion of the Universe 
until their decay rate to quarks is fast enough relative to the expansion rate to allow 



them to evaporate away into baryons and neutralinos [11]. For the case of a non-MAX 
condensate with Q-axitons, the non-linear evolution of a single spherically symmetric 
condensate lump appears to show that /# < 1 once the lump reaches the classically 



stable Q-axiton state [pi ]. It takes some time to reach the Q-axiton state, via emis- 
sion of scalar field waves from a pulsating condensate lump However, in realistic 
models we cannot consider a single, isolated condensate lump, but instead must con- 
sider an ensemble of lumps distributed throughout space. This leads to the possibility 
that scalar waves emitted by one relaxing condensate lump might be absorbed by a 
neighbouring condensate lump. (Since the phase of the Q-axitons (ranging roughly 
from to 29) on average will be the same as the phase 9 of the initial condensate 
lump, absorption and emission of scalar field waves should be relatively efficient and 
without reflection effects due to large (order it) phase differences HSIfl .) As a result it is 
possible that the condensate lumps are maintained in an excited state until expansion 
pulls them far enough apart for absorption/re-emission to be no longer effective. Only 
then can the condensate lumps reach their classically stable Q-axiton state and Jb be 



extracted from lattice simulations |23j of condensate collapse. 

The condition for absorption/re-emission to be important can be estimated as fol- 
lows. We can consider a "lattice" of condensate lumps to form just after condensate 
collapse. These lumps will pulsate about their stable Q-axiton configuration, emitting 
scalar field waves as they do J21] . For simplicity, we will replace the discrete condensate 
lumps with a continuous number density, in order to estimate the number of conden- 
sate lumps encountered by an expanding scalar wavefront. Suppose the condensate 
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collapses to form condensate lumps at t Q , H Q . Consider a scalar field pulse emitted by 
a lump at a later time U. The condensate lump density is 

i 3 

a 



n Q {t) = n Qo (^\ , (7) 



where riq = nQ(t ). The charge absorbed by the surrounding condensate lumps may 
be estimated by considering the area of the outgoing wave removed by encountering 
other lumps. The lumps have area 

Ai » 47rrf , (8) 

where r\ is the radius of the lump. The charge density of the outgoing wavefront of 
radius r is 



Qw , n , 



where Qw is the total charge carried by the wavefront. Thus the charge lost upon 
encountering a condensate lump is ctqAi. The charge absorbed from a wave emitted 
at ti is therefore 

AQ= n Q (t)a Q A l A7rr 2 (t)v(t)dt = n Qi Q w A l v Hr 1 _L ) afda , (10) 

Jti Jai V a J 

where the velocity v(t) of the scalar field waves is initially set by a wavelength of the 
the order of the Q-axiton radius, corresponding to momentum k w 2n\K\ 1 l 2 m [ ^lj| , 
which implies that the scalar field waves are initially mildly relativistic, v(tj) — v « 1, 
and that subsequently t> = v (a,i/a). Therefore the charge absorbed is 



Thus for Oj/a > 1 the charge absorbed is given by 

Qw 5Hi 



(12) 



The condition for total absorption of the outgoing wave is > 1. In terms of if Q and 

/ \ 3/2 

a D , Hi = H Q ( . When the condensate lumps form, the radius of the spherically 
symmetric condensate lumps is approximately the same as the spacing between the 
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lumps, ri « [\K\ l / 2 m) 1 [ 11 , pl|j . Condensate lumps form at H a fa (2\K \m)/a(X ) 
where ct(\ ) ~ 30 [H]. Thus from Eq. ( |T2| ) the condition for absorption by condensate 



lumps to be negligible is that the expansion factor must satisfy a > a c , where 

a, /3^ \ 2 / 3 a 2 / 3 



a 



IK] 1 / 3 ' 



(13) 



This must be satisfied in order to be certain that a numerical simulation of condensate 
fragmentation has allowed the condensate lumps to reach their Q-axiton equilibrium 
configuration, in which case /# can be extracted. In a recent simulation [^3 of conden- 



sate fragmentation on the lattice, Kasuya and Kawasaki found that generally fs = 1 
to a good approximation. Their simulation runs from initial time t = 5 x I0 3 m~ 1 to 
tf = 10 5 m _1 and for values of \K\ = 0.01, 0.05, 0.1. The scale factor at the end of 



their simulation, a/, is thus 

a (l + §# o t) 2/3 ,,,2/3 
•' - • ' 1 ~ ' / » - 6.9 . (14) 



(l + \H t 



2/3 U. 

\ ^7 



On the other hand, the values of a c /a for \K\ = (0.01, 0.05, 0.1) are (31.7, 18.5, 
14.3). Thus the simulation is well within the time scale during which absorption/re- 



emission will be efficient. This might account for the large value of fg obtained in |p23 
as compared with isolated spherically symmetric lump evolution. In order to ascertain 
whether /g = 1 is correct, or whether a \K\ and Q dependent value of fg can occur 



as suggested by single lump evolution pi|, a longer simulation may be necessary. 



The simulation in |23j also shows that for smaller angular velocities of the complex 
condensate field (i.e. amplitude of 02 small compared with <pi), ±Q-ball formation 
becomes important. This is consistent with fluctuations in the phase becoming larger 
than the average phase (02 /<f>i) when the condensate fragments, assuming that the 
initial perturbation of the condensate is unaltered as the angular velocity (02) is re- 
duced. 
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3 Evolution of the Q-axiton 



We next consider if and how a Q-axiton relaxes to the lower energy stable Q-ball. In 
numerical simulations [^1] the Q-axiton is stable on time scales much larger than m -1 ; 
the attractive force between the scalars in the condensate is balanced by the gradient 
energy^] On the other hand, it has been shown that there is no stable spherically 
symmetric breather soliton (such as the Q-axiton) in 3 + 1 dimensions . Thus we 
expect that the Q-axiton will evolve to a lower energy state, most likely a Q-ball. It is 
not apparent from the existing numerical simulation how it loses energy, but from the 
equations of motion we can get some idea of the process. In the equations of motion 
for <pi and cp 2 the argument of the log term, Eq. <p\ + <fi\, is varying periodically in 
time. This is contrast with the case of the Q-ball, for which the amplitudes of 0i and 
02 are the same such that (f>\ + cf)\ is time- independent. It is this property that allows 
one to factor out the time dependence in the Q-ball solution, $(r, t) ~ <p(r)e lujt . For 
the Q-axiton we cannot factor out the time dependence, so we expect the magnitude of 
the Q-axiton field to be varying periodically in time on a time scale m m~ l , albeit very 
slightly based on the numerical results |nj. Thus we expect to find that the Q-axiton 
has a "quivering" solution. This will then result in the emission of scalar field waves 
carrying energy and perhaps charge. The Q-axiton could then either evolve to a Q-ball 
of the same charge or could even, in the opposite limit, completely evaporate away. 
However, since the rate at which energy is being lost by the Q-axiton is very small on 
the dynamical time scale of the scalar field, m _1 , and since the lowest possible energy 
final state would be to have all the charge in the form of a Q-ball, it is likely that 
this will be the final state of the process. If this picture is correct then although the 
Q-axiton appears to be stable on relatively long time scales compared with mT x [ [2T|j , 
on time scales corresponding to the inverse Hubble parameter at the time of Q-ball 
decay it is likely that the Q-axiton will have evolved to the corresponding Q-ball. It 
remains to be shown that Q-axitons do indeed evolve to Q-balls and that the time 



1 Somewhat similar long-lived "pseudo-breather states" are also observed in 3 + 1 Sine-Gordon 
theory, which may represent long-lived pionic states that could be produced in heavy-ion collisions 
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scale for Q-axiton evolution is indeed smaller than that of Q-ball decay. 

The above process is purely classical in nature. It is also possible that the excess 
energy in the Q-axiton could be released via annihilations of the scalars. We can 
consider the Q-axiton to be, to a first approximation, a superposition of two coherently 
oscillating scalar fields and 02- The annihilation rate is n~ 1 drii/dt cx rij (i = 1, 2). 
Therefore since n\ « mcj)\ ^> ri2 ~ m^, the rate at which the <p% amplitude decreases 
will be much larger than that of the (f>2 amplitude. Therefore <pi — > <p2 i- e - the Q- 
axiton tends towards a Q-ball. The final Q-ball charge to the initial Q-axiton charge 
will be Qf/Qi = 4>2f/4>n < (p2i/(pn = Q, the charge of the initial AD condensate. 
Thus if annihilations are able to allow the Q-axiton to evolve to a Q-ball before the 
Q-axiton decays and before the neutralinos freeze out of chemical equilibrium, then 
/b — ► Q fBci, where fsd is the classical value of fs for the Q-axiton. 

We next consider the decay and annihilation of the scalars in the Q-axiton. Q- 
balls and Q-axitons from AD condensate fragmentation have a thick wall profile, well- 
described by a Gaussian profile <p{r) = <ft(0)e~ r2 / R2 , where R ~ (\K\ 1 ^ 2 m)~ 1 is the 
Q-ball/Q-axiton radius [JTTJ, p3fl . Q-ball decay was originally analysed for the case of a 
thin- wall Q-ball with scalars decaying to fermion pairs |H| . The decay rate was found 
to be 

(^] < — (15) 

1 d t . . - 1927T 2 

\ / jermion 

where uj ~ m and A is the area of the Q-ball. This is based on scalar decay filling 
the phase space for fermions of energy m/2 throughout the Q-ball, resulting in Pauli 
blocking of further decay and a Q-ball decay rate which is proportional to A since 
fermions can escape at the surface. In the case of the thick walled Q-ball in the MSSM 
there is no equivalent estimate of the decay rate to fermion pairs. However, we expect 
the order of magnitude to be similar to the thin- wall case for a given 0(0), since the 
decay rate is simply determined by completely filling the fermion phase space within 
the Q-ball volume. In addition, in the case of Q-balls in the MSSM it is possible to 
have decays to pairs of scalars and gauge bosons, which are not suppressed by Pauli 
blocking. In |ll| we considered the decay rate to scalars for the case of a thick walled 
Q-ball. It was shown that most of the Q-ball decay came from scalar decays occuring 



13 



in a region of width Sr pa i? 2 /4r* around r* = (where 7 = Zn 1 / 2 ((70(0) /to) pa 5), 
where tree level decay to pairs of particles coupling to 0, of mass g<j)(r), becomes 
kinematically possible. The resulting Q-ball decay rate was found to be enhanced 
by a factor up to 10 3 compared with the thin-wall fermion decay rate, but is still 
proportional to A. 

We assume that the AD scalar decay rate to pairs of bosons or fermions has the 
form 



a 2 m 3 



r(r) « -T5"> g<p>m 

~ am, g<fi < m , (16) 

(which assumes the largest possible decay rate via exchange of heavy particles once 
g<p > to). We first show that in general Pauli blocking/Bose enhancement effects are 
likely to be significant in Q-ball decay. The maximum number density of particles 
before the occupation number of states becomes greater than 1 is given by 

^p 2 Ap 

Umax ~ 7. — , (17 ) 

where p = |p| = to/2 is the momentum of the final state particles and Ap is the spread 
in momentum of the final state particles. Since the AD scalars in the Q-axiton are 
confined to be within a radius R^ 1 we expect Ap pa R^ 1 and so 

\ K \i/2 m 3 



4tt 2 



(18) 



Pauli blocking/Bose enhancement effects will be important if the density of particles 
from AD scalar decays is greater than n max . Relativistic particles can escape from the 
region of the thick- walled Q-ball 5r around r*, where most of the decays occur, in a 
time r pa Sr pa R/20. Thus the number density of particles from AD scalar decay is 

4nm 3 

Tldecay ~ HI decay T ~ 20|i^| 1 / 2 ' 

where we have used n ~ m0 2 (r) and 2 (r*) ~ m 2 /g 2 . Therefore 

n t \\K\ 
1 h max " I - 1 I 
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This is typically large compared with 1, since \K\ ~ 0.01 — 0.1 in the MSSM [ |I~2]j . 
Therefore Pauli blocking or Bose enhancement effects are likely to be significant in 
the decay of AD scalars in the Q-axiton. The result of this is that all the scalars 
in the region 5r around will rapidly decay. These will be replenished by the Q- 
axiton reconfiguring itself on the dynamical time scale m~ l . In the absence of Bose 
enhanced decay, in a time m~ l the fraction of the scalars in the region 5r which decay 
is Tdecayfn^ 1 ~ oc. Therefore the enhancement in the decay rate is of the order of a" 1 , 
which will increase the decay temperature of the Q-axiton by a -1 / 2 . 

We next compare the decay and annihilation rates. Most of the decays occur 



around r m r* via tree level decays of the AD scalars [H]]. We first show that AD 
scalars in the Q-axiton can annihilate at tree-level to pairs of gauge bosons but cannot 
decay to pairs of gauge bosons. To do this we consider a simple U(l) gauge interaction 
with a pair of "quark" superfields Q\ and Q-i? with U(l) charges ±1. The terms of 
the Lagrangian associated with the gauge interaction are then 

C = D^QlD^Q, + D^Ql^Q^ + (igQ\\^ Ql - igQl^q^ + h.c) 

9 2 12 \n |2\ 2 



+ ^{\QiV-\Q-iV) • (21) 

We can rewrite this in terms of a flat direction scalar $ and its orthogonal partner G, 

$ = -^(g 1 + g_0, (22) 

© = ~7|H2i + Q-i) • (23) 
The kinetic terms in £ are then 
(d^d^ + (0 M e) + 0"e + sfAfA^&Q + 6*6) - gA^i^d^ + iQd^ + h.c.) . (24) 

From this we see that on introducing an expectation value for Re(<&), representing the 
Q-axiton field, the Goldstone boson corresponds to Im(Q). Thus pairs of Q-axiton 
scalars couple to the massive gauge bosons but there is no linear coupling of i?e($) 
to gauge bosons once the Goldstone boson is rotated away. As a result, we can have 
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Q-axiton scalar annihilations to pairs of gauge boson but no decays to pairs of gauge 
bosons. The scalars can, however, decay to gaugino-quark pairs via the interaction 

wQi^q^ - wQ-i^q_^ + h.c. -> -ig$ip e \ + h.c. , (25) 

with decay rate 

r 'decay ~ , (26) 

where a = g 1 jkn. Pairs of AD scalars can annihilate to gauge boson pairs, with 
annihilation cross-section 

a 2 

&ann ~ — K ■ (27) 

Thus the annihilation rate is 

r fl nn = no ann « m(f)(r*) 2 a ann « — , (28) 

where n is the number density of AD scalars. From this we expect that the annihilation 
rate to gauge bosons will be no larger than the decay rate to gaugino-squark pairs. The 
only way to have a significant suppression of the decay rate relative to the annihilation 
rate would be if the decay to gauginos was kinematically suppressed such that the decay 
occurs only to quarks and squarks via small Yukawa couplings. This would depend on 
the mass of the AD scalar relative to the gauginos and the generational structure of 
the mixture of squarks forming the AD scalar. Bose enhancement of the annihilation 
rate to gauge bosons and Pauli blocking of decay rate to gauginos could also help to 
enhance annihilations relative to decays. 

Therefore it is possible that in the limiting case where the non-MAX condensate 
fragments to Q-axitons rather than ±Q-balls, the Q-axiton can evolve to a Q-ball 
before the Q-balls decay, either by classical evolution from the quasi-stable Q-axiton 
state or, less likely, by scalar annihilations dominating over decays. The value of f B 
when the Q-balls decay is particularly small in the annihilation case, which could be 
important in allowing a wide range of neutralino masses to be consistent with the 
baryon to dark matter ratio. However, a more precise analytical and numerical study 
will be required to show exactly how the Q-axiton evolves to a Q-ball. 
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4 Consequences for the Dark Matter and the Baryon 
to Dark Matter Ratio 



4.1 Neutralinos directly from Q-ball decay 



As discussed in the Introduction, a major motivation for B-ball Baryogenesis is the 
similarity, to within an order of magnitude, of the number densities of baryon and dark 
matter particles for the case where dark matter particles have masses of the order of 
mw Q-balls from a d = 6 AD condensate are expected to decay at a temperature 
Td ~ 1 MeV — 1 GeV [11]. Since the scalars in the Q-ball (or Q-axiton) are _R p -odd 



(essentially squarks and sleptons) they will produce one LSP neutralino per scalar. 
For the case of Q-balls from a MAX condensate, assuming that there is no subsequent 
annihilation of the neutralinos from Q-ball decay, the number density of neutralinos 
will be related to the baryon number density by 



or equivalent ly 



n x = 3f B n B 



3/b ('^ ) V B 



Nucleosynthesis constrains fl B ; we consider two possible ranges, 

0.0048 < Q B h 2 < 0.013 



(29) 



(30) 



(31) 



from "reasonable" bounds on primordial element abundances [11 , |34T| and 

0.004 < Q B h 2 < 0.023 (32) 



from conservative bounds |23|, |35||. With h = 0.6 — 0.8 |36|, [LTD we obtain upper bounds 
on the LSP mass neutralino: 



-i (9x\ ( h 

V0.47 U. 



m x < 17.6 ( ^ ) ( — ) GeV 



(33) 



from the reasonable nucleosynthesis bounds and 

? _ 1 f^ x \ ( ^ 



-x<20.8/-^j - GeV 



(34) 



17 



from the conservative bounds. 

For the case of the MSSM, present experimental bounds on the LSP neutralino 
mass depend upon the theoretical assumptions made. The direct experimental lower 
limit on the neutralino mass from ALEPH is m x > 32.3 GeV |P7|| , which requires 
that fs ^ 0.64, assuming the conservative nucleosynthesis limit. If one evolves the 
renormalization group equations assuming universal A-terms and gaugino masses at 
the unification scale then LEP and Tevatron constraints imply that m x > 46 GeV 
(m x > 51 GeV with universal scalar masses also) PU| , which requires that /# < 0.45 
(/b ~ 0.41 with universal scalar masses), assuming the conservative nucleosynthesis 
limits. Thus the case of Q-balls with /b ~ 1 (expected for the MAX condensate), is 
not compatible with MSSM neutralinos. 

If we consider a non-MAX condensate then we can either form ±Q-balls or, as a 
limiting case, Q-axitons. If the Q-axitons cannot evolve to Q-balls before the neu- 
tralinos freeze-out then the extra energy in the Q-axiton relative to a Q-ball of the 
same charge will result in more neutralinos being produced in the decay, such that 
Ib — > Ib/Q in the m x upper bounds. The same correction is necessary in ±Q-ball 
decay. (So as far as the baryon to dark matter ratio is concerned there is no real 
difference between Q-axitons and ±Q-balls in this case.) So for Q < 1 the constraints 
on /b are correspondingly tighter, further disfavouring Q-ball decay as an explaination 
of the baryon to dark matter ratio in the MSSM. (For the ±Q-ball case, /b = 1 is 
expected, whereas for the Q-axiton fg is estimated to be as small as 0.3 for Q « 0.01 
PH| . However, even in this case fs/Q is too large to be compatible with MSSM 



neutralinos.) 

Thus the only way the MSSM can be compatible with direct Q-ball decay to neu- 
tralinos as an explaination of the baryon to dark matter ratio is if we have the limiting 
case where Q-axitons form from a non-MAX condensate and evolve to Q-balls before 
the neutralinos freeze-out. In this case, for small condensate charge Q, /b could be 
small enough to allow dark matter MSSM neutralinos compatible with experiment to 
come directly from Q-ball decay, especially in the case where annihilations are respon- 
sible for the evolution of the Q-axitons to Q-balls. 
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If we take the view that the baryon to dark matter ratio is indeed due to decaying 
Q-balls then, modulo the possibility of Q-axiton formation and evolution, we must 
consider an alternative to MSSM LSP neutralinos. A natural possibility is to consider 



the NMSSM with a mostly singlino LSP gf| |37|, g|. With f B « 1 and Q < 1, Q- 
ball decay imples that the LSP must have a mass less than about 20 GeV. Current 
experimental bounds combined with renormalization group evolution indicate a lower 
bound of 3-5 GeV on the singlino LSP mass [p8|| . (Singlinos may be excluded by 
vacuum stability considerations, but this depends on the universality conditions at 
the unification scale [P^Q. 

So, given that J'b/Q ^ 1 is likely when the Q-balls decay, Q-ball decay combined 
with the baryon to dark matter number ratio and experimental bounds on neutralino 
masses rules out the MSSM but is compatible with the NMSSM with a mostly singlino 
LSP of mass less than about 20 GeV. This may be regarded as a prediction of the 
B-ball Baryogenesis explaination of the baryon to dark matter ratio. 

4.2 Effect of annihilations 

In the above it has been assumed that the neutralinos from Q-ball decay do not 
annihilate. However, if a large enough density of neutralinos is produced in Q-ball 
decay, they will annihilate with each other. In this case the number density of LSP 
neutralinos is determined by the annihilation cross-section and is not directly related to 
the baryon number density. It will also be spatially constant, ruling out the possibility 
of neutralino isocurvature fluctuations correlated with the baryon number |38| . 



The upper limit on the number of LSPs at a given temperature is [TT 



n Umi t{T) = ( — ) , (35) 

V< CTV >annJ T 

where a is the annihilation cross-section, v is the relative velocity of the neutralinos 
and < ... > denotes the thermal average. For neutralinos < av > an n= a + bT/m x , 



where a and b are determined by the parameters of the SUSY model ||39|| . For the 



case of light neutralinos m x < m\y, the a term is negligible and the annihilation 
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cross-section is b dominated, in which case nu m i t oc g(T) l l 2 T. Therefore 

(9(Tf)Y /2 (T f \\ (T) 
[~^Tj) [t) UreUc{Th 

where 



1/2 

nu m u(T) = I ^J^T ) ( 7F ) n relic(T) , (36) 



g(T)\[T\ 3 / H 



n "' UT) = IjMJ (37) 

is the thermal relic neutralino density and Tf is the neutralino freeze-out temperature 
(typically Tf « m x /20). Thus in order to account for the baryon to dark matter ratio 
the LSP density from Q-ball decay must satisfy 

n re iic{Td) < n LSP (T d ) < n Hmit (T d ) . (38) 

In general 

n LSP (T d ) = . (39) 



Thus the condition riLSp(Td) < nu m it{T d ) implies that |TT 



MJ V/J \m) m (ioo^v)' (40) 



Since the decay temperature is estimated to be between 1 MeV and 1 GeV ||11|| , for 
reasonable values of Q and fs (say Q ^ 0.1) we expect that it is quite natural for 
this to be satisfied, especially for the case of a maximally charged condensate with 
Q = f B = 1. However, if T d > T dc then the LSP density is given by n Umit {T d ), 
independently of the details of Q-ball decay, 

(T \ ^ 
-jrj n L sp(T d ) . (41) 

Although the direct connection between the baryon number and dark matter number 
density is lost, the neutralino dark matter will still be non-thermal in nature. 



5 Non-Spherical Collapse and Velocity Effects 

In the above discussion we have been considering values of /# based on the spherically 
symmetric collapse of a single condensate lump and the evolution of an ensemble of 
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static lumps being pulled apart from each other by the expansion of the Universe. 
However, it has been suggested that the effects of the velocity of condensate lumps 



after fragmentation [40[ and the non-spherical collapse of condensate lumps |23| , [40f 
may have very significant effects on the evolution of the ensemble of condensate lumps. 
In a recent numerical simulation |^D| it was observed that after condensate fragmenta- 
tion the lumps initially have a mildly relativistic velocity and can undergo collisions. 
(However, the lump velocity may depend upon the assumed form of the initial per- 
turbation at H pd m. For example, a standing wave perturbation may result in static 
lumps.) In addition, the large lumps which initially formed subsequently fragmented 
into ±Q-balls. One question here is whether the formation of ±Q-balls is due to the 
dynamics of condensate collapse or due to the linear growth of the perturbations of 
the charge density, as discussed in Section 2? As well as mildly relativistic large Q- 
balls, smaller relativistic Q-balls were formed during condensate collapse. These, it is 
suggested, can effectively thermalize the distribution, allowing an analytical treatment 
[fl(J]. The possibility of forming a significant number of small Q-balls in condensate 
fragmentation could have a significant effect on fs- These could be emitted during 
the evolution of a non-spherical condensate lump, either due to being initially non- 
spherical after formation or due to the fusion of two colliding condensate lumps; for 



example, in |23[ it was shown that in a 1 + 1 dimensional simulation of Q-ball fusion, 
7% of the charge was emitted in the form of very small Q-balls (or possibly scalar 
field waves, as the lattice simulation cannot distinguish these cases [fffil). Since the 
rate of decay of Q-balls is proportional to Q _1 and the decay temperature oc Q -1 / 2 , 
it is possible that very small Q-balls could decay before the neutralinos freeze out of 
chemical equilibrium, Td > Tf. In this case, since the neutralinos from Q-ball decay 
thermalize to become part of the thermal relic density and so are not related to the 
baryon number from the decaying Q-balls, the charge in small Q-balls is effectively 
removed from the Q-ball ensemble before the larger Q-balls decay, reducing fs ■ So the 
tendency of non-spherical lumps and velocity on the formation and evolution of con- 
densate lumps is to reduce J'b as compared with the idealized spherically symmetric 
static case. 
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6 Conclusions 



We have considered a number of issues connected with the fragmentation of the Affleck- 
Dine condensate and the formation of late-decaying Q-balls, in particular the fraction 
of the baryon number in the Q-balls when they decay, For a non- maximally charged 
condensate we expect that either ±Q-balls or Q-axitons will form, depending on the 
initial fluctuation of the phase of the condensate at H m m. We showed that in the 
limiting case where Q-axitons form numerical simulations of condensate fragmentation 
must be run for a sufficiently long time to ensure that absorption/re-emission of scalar 
field waves does not maintain the condensate lumps far from their classically stable 
Q-axiton state, so allowing the value of fs to be extracted. 

A major motivation for B-ball baryogenesis is a natural explaination for the baryon 
to dark matter number density ratio. In general, Q-ball formation with f B = l 
from fragmentation of a MAX condensate and ±Q-ball formation from a non-MAX 
condensate produces too many neutralinos per baryon number to be consistent with 
experimental bounds on MSSM neutralino masses if we wish to account for the baryon 
to dark matter particle number ratio via Q-ball decay. The only possibility is to have 
the limiting case of Q-axiton formation from a non-MAX condensate, with evolution 
of the Q-axitons to Q-balls before neutralinos freeze out of chemical equilibrium, such 
that fg 1 is possible. This is particularly possible if the Q-axiton evolves to a 
Q-ball via annihilations. 

It is also possible that more complex features of AD condensate fragmentation 
which can only be studied numerically, such as the formation and evolution of non- 
spherically symmetric condensate lumps and collisons of moving lumps, might reduce 
the value of /#, but it remains to be seen if these effects can reduce it sufficiently to 
accomodate MSSM neutralinos in the typical case of Q-ball formation from a MAX 
or ±Q-ball formation from a non-MAX condensate. 

In the typical case with fs/Q ^ 1 (rather than the limiting case of Q-axiton 
formation followed by evolution to Q-balls), the LSP neutralino mass from direct Q-ball 
decay must be less than about 20GeV. This is too large to be consistent with MSSM 
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neutralinos (m x > 32.2 GeV experimentally), but is consistent with singlino dark 
matter in the NMSSM. Thus the B-ball baryogenesis explaination of the baryon to dark 
matter ratio strongly suggests that dark matter is non-thermal and that the LSP should 
be associated with an extension of the MSSM with an experimentally allowed light 
LSP, the NMSSM being the most obvious possibility. In addition, depending on details 
of the inflation model, there may be observable isocurvature density fluctuations due 
to the transfer of baryon isocurvature fluctuations to neutralinos via Q-ball decay [ |31fj . 
Although all the evidence for B-ball baryogenesis is circumstantial, and each element 
could have an alternative explaination, if it turns out that the NMSSM is realized 
in nature with a light singlino LSP, if the LSP dark matter density is inconsistent 
with a thermal relic density and if isocurvature density perturbations are observed, 
then d = 6 B-ball baryogenesis would become a favoured scenario for the simultaneous 
origin of both the baryon number and dark matter density. 

Finally, it should be emphasized that much work remains to be done to clarify 
the physics of unstable AD condensates and B-ball Baryogenesis. Given that d = 6 
AD baryogenesis is a serious candidate for the origin of the baryon asymmetry in 
SUSY models and that late-decaying Q-balls could naturally account for the baryon 
to dark matter ratio, it is important to develop in detail the physics of AD condensate 
fragmentation in the early Universe. 
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